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Abstract 

We give a construction for a self-test for any connected graph state. 
In other words, for each connected graph state we give a set of non- 
local correlations that can only be achieved (quantumly) by that par- 
ticular graph state and certain local measurements. The number of 
correlations considered is small, being linear in the number of vertices 
in the graph. We also prove robustness for the test. 

1 Introduction 

Self-testing is a process where a skeptical classical user attempts to verify the 
operation of a collection of quantum devices without trusting any of them a 
priori. Importantly, we wish to make as few assumptions as possible about 
the operation of the devices and in particular we do not bound the dimen- 
sion of the state space for each device. However we do make the necessary 
assumption that the quantum devices are not allowed to communicate with 
each other. Despite these severe restrictions on our knowledge it is possible 
to devise self-tests for a number of different situations. 

Self-testing was first introduced by Mayers and Yao |MY04] who described 
a self-test for a maximally entangled pair of qubits (EPR pair) along with a 
small set of local measurements. Meanwhile, self-testing of gates was intro- 
duced by van Dam et al. |vMMS00] in the scenario of known Hilbert space 
dimensions. These two results were extended to testing of circuits over a real 
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Hilbert space by Magniez et al [MMM0Q6]. Most recently, McKague and 
Mosca MM 1 T[j reproved the Mayers- Yao result and extended it to allow for 
testing of a larger set of measurements including measurements over the full 
complex Hilbert space. 

In this paper we use proof techniques developed in [MMlOj to define self- 
tests for the graph state for any connected graph. This family of self-tests 
is efficient in the number of measurement settings, requiring only two or 
three measurement settings on each vertex, depending on the graph. As well 
the total number of correlations tested is small, only one per vertex plus an 
additional 3 at most. We also prove that the self-tests are robust. 

1.1 Graph states and notation 

A graph G is composed of two sets: a set V of vertices, and a set E C V x V 
of edges. For our purposes we suppose that (v,v) E and (v,u) G E 
whenever (u, v) G E. Two vertices u, v are said to be adjacent if (u, v) G E. 
A cycle is a sequence of vertices in which each vertex occurs at most once, 
each vertex in the sequence is adjacent to the next vertex in the sequence, 
and the last vertex is adjacent to the first. A subgraph G' of G is a graph 
(E', V) with E' C E, V C V. An induced subgraph is a subgraph in which 
E' = {(u, v) G E\u,v G V'}, so the subgraph contains all edges between 
vertices of V in the original graph. The neighbours N v of a vertex v are the 
vertices to which v is connected with an edge, i.e. N v = {u G V\(u, v) G E}. 
A bipartite graph is a graph in which the set of vertices may be partitioned 
into two sets S and T, each of which has no edges within it. So the induced 
subgraphs on S and T have no edges. An important property of bipartite 
graphs is that they are exactly the graphs which contain no cycles with an 
odd number of vertices. A graph is connected if for each pair of vertices u, v 
there is a sequence of adjacent vertices beginning with u and ending in v. 
For more detail regarding graph theory see Diestel |DielO] . 

A graph state consists of a set of qubits indexed by the set of vertices V, 
each prepared in the state \+) v = ((0)^ + |1) J, followed by (CTRL — Z) uv 
operations between pairs of qubits where the corresponding vertices u, v in 
the graph are adjacent. If the graph is not connected then the graph state 
will be a product state of graph states on the separate components. Hence 
connected graphs form the interesting cases. 

Graph states are also characterized by their stabilizer group. Let the 
operators X v and Z v be the Pauli operators X and Z applied to qubit v, 
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tensor product with I on all other qubits. If P is a Pauli and S C V then 

p s =u p - 

The stabilizer group for a graph state on the graph G = (V, E) is generated 
by 

S v = {X v Z N «\veV}. (2) 

That is, for each vertex v there is a stabilizer operator with X operating 
on v and Z operating on each of i>'s neighbours. Note that there are n 
such operators, they pairwise commute and are independent. Hence there is 
exactly one state with this stabilizer group. That is to say, the graph state 

is the unique state for which S v = for each v G V. 

As one additional piece of notation, we will frequently need to deal with 
products of stabilizers on a subset of vertices. For this case we define 

Z N(S) = "Q Z N V (3) 

where the factor Z v appears if v has an odd number of neighbours in S. 
1.2 Self-testing definitions 

Consider the following black-box scenario: we are given a set of devices, each 
with a knob labeled with a number of settings, a pair of lights labeled ±1, 
and a button. After we select a setting and push the button one of the 
lights turns on. We are told that the devices jointly share a state which is 
measured, according to the knob setting, in a specified basis. Our goal is 
to determine if the black-boxes are operating according to their specification 
using only the external controls of the boxes. Additionally we may isolate 
the boxes to ensure that they do not communicate. 

We begin with a reference experiment consisting of an n-partite system in 
the state together with local measurement observables M^ m on subsystem 
j with measurement setting m G {0,1, ... ,kj}. The measurement setting 
m = corresponds to no measurement, which we may represent with the 
identity. The reference experiment represents the specification for how the 
black-boxes supposedly operate. In particular, we assume that the state and 
observables are known. 
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In addition, we have a physical experiment consisting of an n-partite phys- 
ical system in the stated \ip') together with local measurement observables 
M'j m on subsystem j, with m G {0, 1, . . . , kj}. Again we may take Mj = / 
indicating that we do not measure the subsystem. We place no bound on 
the dimension of the Hilbert space of each subsystem, but assume that it 
is finite. The physical experiment represents how the black-boxes actually 
operate. 

If a physical and reference experiment have the same number of subsys- 
tems and the same number of measurements on each subsystem, then we 
say that they are compatible. Note that we will always deal with the case of 
two-outcome measurements, so that all observables have eigenvalues ±1. In 
principle, though, the definitions can be extended to other types of measure- 
ments. 

To be more specific about our task, we introduce two notions, simulation 
and equivalence. 

Definition 1. Let a physical experiment and a compatible reference exper- 
iment be given as above. We say that the physical experiment simulates 
the reference experiment if for each measurement setting m = (mi, . . . , m n ), 
nij G {0, . . . , kj} we have 

n n 

W\ (8) Km, W) = (4>\ ® M j>mj |V) • (4) 
3=1 3=1 

Here it will be sufficient to consider only a subset of possible measurement 
settings. In this case we include the measurement settings of interest in our 
description of the reference experiment. 

Definition 2. Let a physical experiment and a compatible reference experi- 
ment be given as above. We say that the physical experiment is equivalent 
to the reference experiment if there exists a local isometry 

$ = $i®---®$„ (5) 

and a state \junk) such that, for each j , and m G {1, . . . ,kj} 

= \junk)®\i>) (6) 
$(il4 m |V0) = \junk) ® M hm \^) (7) 

1 We consider only pure states, but since the Hilbert space of the physical system has 
unbounded dimension we may easily add a purification to mixed states. 
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where \junk) is in the same Hilbert space as \ip') . 

When describing any physical system we must first fix a reference frame, 
and decide which components to describe and which to leave out. Thus we 
may take a description and apply local changes of basis, or add ancillas and 
arrive at another, perfectly acceptable, description of the system. These 
two operations are invisible from the perspective of classical interactions 
with devices so we can never rule them out. This motivates our definition 
of equivalence, which takes such ambiguities in quantum descriptions into 
account. 

Throughout the remainder of this paper we will used primed d^'), X', 
S' v etc.) to denote physical measurements and states and unprimed for ref- 
erence measurements and states. Note that S' v = X' v <g> Z' N ^ and other 
derived physical measurements are defined in terms of the local physical 
measurements. Also, although we use the letters X and Z for the physical 
measurements, these need not be Pauli matrices, and we assume nothing 
about their structure other than what we mention explicitly. 

1.3 Main results 

A self-testing theorem specifies a particular reference experiment and states 
that if a physical experiment simulates the reference experiment, then it is 
equivalent to it. That is to say, for a particular experiment simulation implies 
equivalence. Our main result is to show that this is the case for the following 
two reference experiments. 

Definition 3 (Reference experiment 1: connected graph with an odd induced 
cycle). Let G = (V,E) be a connected graph containing an odd induced cycle 
C = (V',E'). Let be the corresponding graph state with stabilizers S v , 
v £ V. The reference experiment consists of the state the stabilizer 
measurements S v and the measurement X v Z N( - V ) . 

It is easy to show that a graph which contains any odd cycle contains 
an induced cycle. Thus reference experiment 1 is applicable to all connected 
non-bipartite graphs. 

Definition 4 (Reference experiment 2: connected graph). Let G = (V,E) 
be a connected graph with at least two vertices. Let \tp) be the corresponding 
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graph state with stabilizers S v ,v G V. Choose a fixed edge (u,v) G E and 
define 

D u = ^ (X u + Z u ) (8) 

The reference experiment consists of the state \tp) , the stabilizer measure- 
ments S v and the measurements 

Z' U Z ,N " (9) 

D U Z N - (10) 
D u X v Z N ^ {u] (11) 

In appendix [A] we show that the D measurements are required since for 
a bipartite graph all measurements using X and Z alone can be simulated 
using a classical hidden variable model. 

Theorem 1. If a physical experiment is compatible with reference experiment 
1 (2), and simulates it, then the physical experiment is equivalent to reference 
experiment 1(2). 



2 Proof of main result 

The proof consists of three sections. First we determine the expected values 
for the measurements in the reference experiment. Next we show that if the 
physical experiment simulates the reference experiment then the X' and Z' 
operators anti-commute. Finally we construct the local isometry and use the 
anti-commuting property of the X' and Z' operators to show equivalence. 

2.1 Probability distribution from graph states 

We first derive the probability distributions that arise from a graph state 
with trusted measurements. This establishes the conditions that a physical 
experiment must meet in order to simulate the reference experiment. 
Clearly, the stabilizer measurements all satisfy 

{iJ>\S v \i/>) = l. (12) 
For reference experiment 1, we need one additional measurement. 
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Lemma 1. Let G = (V, E) be a graph and let \ip) be the corresponding graph 
state. Let V C V and let G' = (V',E') be the induced subgraph on V . If 
each v € V has even degree then 

(-1)\ E '\X V 'Z N ^ \i/>) = |V> (13) 

Proof. Consider the product 



(14) 



First note that there will be an X v factor for each v 6 V . As well, there will 
be a Z u factor for each v G V adjacent to u. Canceling pairs we see that 
there will be an overall Z u factor exactly when there are an odd number of 
neighbours of u in V. Hence the Z factor will be Z N ( V >. We only need to 
determine the sign. Note that the Z u , u ^ V factor all commute so we need 
not consider them any more. 

The order of multiplication in equation (fT4"j) does not matter since the 
stabilizers all commute. For convenience, then, we may write the product 
as the product of the rows of a matrix with each column corresponding to a 
v £ V and each row a stabilizer. We choose the order of the rows so that 
the Xs appear along the diagonal For a 5-cycle, for instance, we have 

X Z I I z 
Z X Z I I 

I Z X Z I. (15) 
I I Z X Z 
Z I I Z X 

The factor on each vertex equals the product of the entries in the correspond- 
ing column. In each column there is one X and one Z for each neighbour. 
The factor will be either ±XZ or ±X, depending on whether there is an 
odd or even number of Zs. The sign depends on the number of Zs above 
the X, since we must use the fact that XZ = —ZX once for each such Z. 
Combining the signs from all vertices, there is a —1 factor for each Z above 
the diagonal, and hence one for each edge in G'. The overall sign, then, is 

(-i)in □ 



2 The matrix may be constructed by taking the adjacency matrix of G', which has a 1 
in the u,v position when (u,v) E E' , replacing the diagonal with Xs, the Os with Is and 
the Is with Z. 
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For reference experiment 1 we consider an odd cycle, and hence we obtain 

(ip\x v 'z N( - v y\i/j) = -i. (i6) 

Reference experiment 2 has three measurements other than the stabilizer. 
First we have Z u Z Nu . This is just S u with X u replaced by Z u . Since X and 
Z anti-commute we have 

M Z U Z N « |V) = 0. (17) 
From this, and linearity, we obtain 

WD U Z»-\1,) = ±=. (18) 

Finally, the operator D u X v Z Nv ^ u ^ is a linear combination of S v and S v with 
Z u replaced with X u . As above, then, we find 

<V|AA^M|^ = _1 (19) 

2.2 Statistics imply anti-commuting observables 

We now suppose that the physical experiment simulates either reference ex- 
periment 1 or 2 and show that this implies that the X' and Z' measurements 
on each vertex anti-commute (on the support of |^)). 

First, note that (if)'\ S' v \ip') = 1 implies S' v \ip') = and similarly for 
other measurements. This allows us to immediately drop probabilities and 
deal with states directly. 

As a first step towards our goal, we prove a type of induction lemma 
which says that if the X' and Z' observables anti-commute on vertex, then the 
same is true for an adjacent vertex. Thus we need only show anti-commuting 
observables on one vertex, and apply the lemma repeatedly along paths to 
all other vertices (since G is connected.) 

Lemma 2. Given a graph G with (u,v) G E. If observables X f v , Z f v ,X' u , Z' u , 
and {Z' w \w G N u U N v } and state \ip') satisfy 

s' u W) = s' v W) = W) (20) 

{X'Z') v \^) = -{Z'X') v \i,') (21) 

thcfi 

(X'Z') u m = -(Z'X') u m (22) 



8 



Proof. From the fact that (u, v) G E we obtain 



IVO = (M%^5XW> (23) 

= {Z'X') U X' U Z' V X' V Z' U X' U Z' V X' V Z' U W) (24) 

= {X'Z'UZ'X'UZ'X 1 )^ 1 ) (25) 

= -(X'Z')„(Z'X') B (X'Z') W |^'> (26) 

= -{X>Z') u \i>>) (27) 



(28) 
□ 

For reference experiment 1 we show that the observables X' and Z' anti- 
commute for each vertex in the induced odd cycle. 

Lemma 3. Let G = (E, V) be a connected graph and let C = (£", V) be 
an induced odd cycle of G and let u G V . If observables X' u , Z' u for u G V , 
{Z' w \w has a neighbour in C} and state \ip') satisfy 

S' u W) = W) (29) 

- X' V 'Z' N ^ = W) (30) 
Then (X'Z') U = -(Z'X') U for each u G V . 

Proof. Number the vertices in the cycle 1 through k so 1 is adjacent to 2, 
etc.. Without loss of generality we may assume that u is vertex 1. We next 
consider the following state: 

- X ' v 'z' N ^ n s' 2j n ^,--1 w) = w) (si) 

3=1 3=1 

Note that the factor Z' N(yVn> is cancelled by Z operations arising from the 
products of the S' v . We may write the product as the product of the rows of 
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the following matrix, where column j corresponds to vertex j in the cycle: 



-X' 
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X' 


X' 
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X' 


X' 


Z' 


X' 


Z' 
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.. I 
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Z' 


X' 


Z' . 
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I . 


.. Z' 


X' 


Z' 


X' 


Z' 
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I . 


.. I 
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Z' 
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Z' 


X' 


Z' 


I . 


.. I 
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I 


Z' 


X' . 


.. I 


I 


I 


Z' 


I 


I 


I 


I . 


.. I 


Z' 


X' 



(32) 



In each column there are two X' operators and two Z' operators. Also, 
their arrangement is such that, for every column except the first, the two X' 
operators are next to one another, so they cancel directly, and similarly for 
the Z' operators. Hence 

- (X'Z') U (X'Z') U = W) (33) 
The desired result follows immediately. □ 

For reference experiment 2, we have one additional measurement on a 
particular vertex u. We use this extra measurement to establish that the X' 
and Z' measurements on u anti-commute. 

Lemma 4. Let G = (V,E) be a connected graph with (u,v) G E. If observ- 
ables D' u , X' v , Z' v , X' u , Z' u , {Z' w \w G N u U N v } and state \ip') satisfy 

S' u W) = S> W) = W) (34) 



m z' u z' n - w) 

W\Df u Z' N -W) 
{^\D' U X V Z' N ^ W) 





1 
1 



(35) 
(36) 

(37) 
(38) 



then-{X'Z%W) = {Z'X') u \i,') 
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Proof. Since X' U Z' N - = 1 we have X' u = Z' n "\^j). Similarly, 
Z' u \if>') = X' V Z' N ^ U \^'). Along with (if>'\ Z' U Z' N " \^') = we find that 
X' u \ip') is orthogonal to Z'Jip'). We also obtain D' U Z' U \ip'} = ^= and 
(ip'\ D' U X' U |^') = ^. Since D(, |^') has norm 1, we find 

D' U W) = ^=X' U W) + Z' U ^=W) (39) 

Further, since {D' u f = I = (Z'J 2 = (X' u ) 2 , and 

W) = {D'ufW) (40) 

= -^D' u (Z' N -+X v Z' N v\ u )W) (41) 
v2 

= 1 -{Z' N -+X V Z' N ^){X' U + Z' U )W) (42) 

= i (2/ + (X'Z% + (Z'X%) W) (43) 

(44) 

In order for this to be true, we must have 

(X'Z') U \^) = -(Z'X') U \^'). (45) 

□ 

We conclude with a technical lemma that allows us to exchange X' v op- 
erations for Z' v operations. 

Lemma 5. Let G = (V,E) be a connected graph and let X' V ,Z' V for v G V 
and (and D u for some u G V) be a physical experiment that simulates 
reference test 1 (or 2). Let G' = (V',E') be an induced subgraph of G. Then 

\$) = Z' N ^ W) (46) 

Proof. We use the previous lemmas to conclude that X' V Z' V \ tp') = —Z' V X' V \ip') 
for each v. Then we repeat the argument used in the proof of lemma [U 
Essentially, we look at the product 
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Writing this product out as a the product of rows of a symmetric matrix with 
X's along the diagonal, we see that in order to get all the X's together we 
must use the anti-commuting relation once for each Z' above the diagonal. 
Since there is one Z' above the diagonal for each edge, we obtain the factor 

(-i)in 

□ 



2.3 Constructing the isometry 

The local isometry $ that we use to show equivalence between the physical 
experiment and the reference experiment is the tensor product of isometries 
<& v for various v G V, is in the circuit shown in figured) 
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\input) v 



H 



H 



Z' 



XL 



Figure 1: Circuit for 

The circuit is based on the argument used by Mayers and Yao in their 
original EPR test. It may be seen as a type of SWAP gate, decomposed 
into three CNOT gates. Here the first CNOT gate is omitted since the 
ancilla is always initialized in the state |0). The Hadamards and Controlled 
Z operation replace a CNOT targeted on the ancilla. With these points in 
mind, we see that when Z' v and X' v are indeed qubit Pauli operators the 
circuit defines a SWAP operation. 

We will now calculate the result of <3> applied to 

aKIVO) = ^E<SW (/ + (-1)""^) W) \x) (48) 

with x = (x v ) ve y G {0, l}l y L Applying the anti-commutation relation, this 
simplifies to 

$ (^» = E <8> ( J + Z 'v) X '* V l^> l*> • ( 49 ) 

x uev 

Using lemma [5] and the fact that (/ + Z' V )Z' V = I + Z' v we finally find 



uev 




^E(- 1 ) e(I ' ) i-)) 



(50) 
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where e(x) is the number of edges in the induced subgraph on the set V x = 
{v G V\x v = 1}. 

Set |0) = ^7= J2 X (~ l) 6 *^ \x). Consider S v \x) for some x. This will be 
± \x © 1„) where 1„ is the binary vector with 1 in position v and everywhere 
else. The sign may be computed as follows: for each Z u component of S v , if 
x u = 1 a —1 factor will be introduced. This happens when (u,v) G E and 
u is in V x . We may see this as either removing or adding the vertex v and 
adding a —1 factor for each edge between v and another vertex in V x . Thus 
S v (— l) e w | a;) = (— i)^© 1 ") \ x © l v ). In other words, this exactly produces 
the correct sign on each \x) so that S v \<f>) = \<p) and in fact \<p) — 

Now consider &(X' V |V')) for some v. After anti-commuting the X' oper- 
ations we have 

hk m = l £ (g) (/ + x:-x; ivo I*) . (si) 

In this equation, we may simply replace X' v Xv X' M with X' v Xv ® lu , where l u is 
the vector with Os everywhere, except position u. After applying lemma [5] 
we arrive at 

HK W)) = ( i (/ + z;) |v') ) ^(-i)^©i«) \ X ) . (52) 

A change of variable, x h- >■ x © l u , and the fact that X u \x) — \x © 1 M ) gives 
the final result, 



$(A: |V')) = ^= <jg) (I + 4) |V') j X„ |V>) . (53) 
A similar analysis shows that 

W IVO) = ® ( J + 4) iVol ^ IV>) • (54) 

Recall from the proof of lemma H] that D' v may be written as D' v \ip') = 
i ( A; + Z£) |V')- By linearity, then 

HD> W)) = (-L 0(1 + Z>) |V') ] D v |V) • (55) 
This concludes the proof of theorem [TJ 
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3 Robustness 

In this section we will show that the main theorems are both robust. 



3.1 Definitions and main theorem 

First, we modify the definitions of simulation and equivalence to allow for 
small deviations from the reference experiment 

Definition 5. Let a physical experiment and a compatible reference exper- 
iment be given as above. We say that the physical experiment e-simulates 
the reference experiment if for each measurement setting m = (mi, . . . ,m n ), 
rrij G {0, . . . , kj} we have 



3=1 3=1 



3<"H 



< e. 



(56) 



Definition 6. Let a physical experiment and a compatible reference experi- 
ment be given as above. We say that the physical experiment is 5-equivalent 
to the reference experiment if there exists a local isometry 



and a state \junk) such that, for each j , and m e {1, . . . ,kj} 



\®(W)) - \junk) <g> 



< 5 



(57) 



(58) 
(59) 



where 5 = 15n2 + 5n ^J~l ($ —9J and \ junk) is in the same Hilbert space as ■ 

Theorem 2. Let a graph G be given with \V\ — n. If a compatible physical 
experiment e-simulates reference experiment 1 (2) then it is 5-equivalent to 
it with b=\ (5n 2 + \\n + 4) yji (5 = (2n 3 + An 2 + n)y/e + 13(|n 2 + n)ei ) 



3.2 Proof for reference experiment 1 

First we note that if M > 1 — e then 



(60) 
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Next, suppose that we have \\\^) - M < a and \\\tp) - N \tp)\\ 1 < (3. 

Using the triangle inequality and the fact that ||M|| = 1 we have 

\\\$)-MN\il>)\\i<<x + P- (61) 

The remainder of the proof will use these estimations repeatedly, along with 
the triangle inequality. We need only count the number of times this happens, 
which is the same as the number of operators multiplied together. 

First, for lemma [3] let c be the size of the induced cycle. We multiply 
c + 1 operators together. Thus we conclude that for a vertex u in the induced 
cycle 

\\(X'Z% |V/> + (Z'X% 1^)11! < 2(c+ l)vfc (62) 

Next, for lemma|2]we multiply four operators, then invoke the anti-commuting 
property on one of the vertices. This gives 

1 1 (X'Z') U W) + (Z'X% IV) || x < 8v^ + /3 (63) 

where (3 is \\{X' Z') v \ip') + (Z'X') V \ip')\\ v v being neighbouring vertex. We 
may apply lemma |2] along paths from vertices in the induced cycle in G. Let 
I be the length (number of edges) of the longest path. Then for any vertex 
u we find, at worst, 

\\{X'Z') U + (Z'X') U |^)Hi < 2(4/ + c + l)V~e. (64) 

Lastly, for lemma |5l we multiply \V'\ operators, and apply the anti- 
commuting relation \E'\ times. Thus 

-1)^X' V ' - Z' N ^ |V') 1 1 < 2 + (4/ + c + 1)\E'\) yft. (65) 

We are now ready to analyze the proof of the main theorem for reference 
experiment 1. To arrive at equation H9] we apply the anti-commutation rela- 
tion. This happens once for each 1 appearing in x, for each possible x, for a 
total of n2 n ~ l times. We may find this by pairing values x and x © 111 ... 1. 
There are 2 n_1 such pairs and each pair contains n Is all together. Multiply- 
ing by the normalization factor we find 



$ (i^'))-^E(8)( / +^) x ^ i^>i*> 



x vev 



< n(4/ + c+ l)Ve- (66) 
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For equation [SU] we use lemma \5[ once for each possible value of x. Again, 
the estimate depends on the number of Is in x, summed over all possible xs. 
As well, it depends on the number of edges in the induced subgraph. An 
edge (it, v) will be counted only when x u = x v = 1. This occurs for 1/4 of all 
xs. Summed over all possible xs and edges, then, the number of times edges 
are counted is 2 n ~ 2 \E\. Again, we multiply by the normalization factor 
This gives our final estimate: 

V tiev / x 

< (n(4Z + c + 1) + {n + (4/ + c + l)|£|/2)) (68) 

= K4Z + c + l)(n + ^)+nJ (69) 

where e(x) is the number of edges in the induced subgraph on the set V x = 
{v G V\x v = 1}. 

Note that when calculating $ (X' u \ip')) etc. we did not use any more 
estimations, we simply rearrange when lemma [5] is applied. Thus the same 
robustness applies. 

As a last estimation, we note that I and c cannot be larger than n, and 
\E\ < n 2 . We may thus set 8 — | (5n 2 + lira + 4) i/e. 

Note that we may make much better estimates if some properties of the 
graph are known. For example, if every vertex lies in a triangle and the 
max degree is 6, as in the case of a lattice of triangles, we may instead set 
5 = 17n-y/e. 



(67) 



3.3 Proof for reference experiment 2 



Much of the same analysis may be used for experiment 2. Indeed, since 
the only difference in the proofs for the non-robust results is how the anti- 
commuting property is proved, we may simply replace the estimation for 
lemma |3] with that of lemma HJ 

We begin, then, with e-simulation and prove a robust version of lemma SJ 

xt. +Z' i _ i 



First we wish to estimate a = D' u \ip : 
from the previous section, we have 

\\KW)-z' N "W\ 
\\z' u W)-x' v z' N ^W)\ 



V2 



< 
< 



Using techniques 



2^1 (70) 
2^. (71) 
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These along with the triangle inequality give an upper bound for a of 

Z'N U _|_ x' Z' Nv \ u 



2V2~e + 



Dl. 



V2 



(72) 



Expanding the second term, we get 



Z'Nu _|_ x' z iN v\ u 



' 1 + = '-^ | V) - V2 ((VI D> U Z' N « |V) + ^D^Z^A" |V)). 

(73) 

Since | \Z' U |V) - X' V Z' N -\ U |V)| | x < 2^ and ||Z ,JV « (V)^ = 1 we find 

I (V| Z IN -Z' U |V) - Z /7V "X^ /A ^ U |V) I < 2v^- (74) 



By hypothesis, 



Z' N "Z' u \r>\ <e, so 

Z' N u+X'Z ,N v\v 



(VI Z ,JV "X;Z ,JV "\ U |V) | < 2^ + e. 
2 = 1 + Re (VI Z' N -X' V Z' N ^ U |V), 



Meanwhile /3 2 = |V) 

so |1 -/3 2 | < 2^1+6. 

Finally, by hypothesis |(V| I^Z /JV « |V) + D' u X' v Z lN ^ u |V) - < 

2e. Combining these facts we find a < 2^ + \j2y/e + (1 + 2<\/2)e. 
Now we wish to estimate 



( X 'u + Z 'u. 



By the fact = 1 we have {D' u f |V) - K^pt W) 

larly, since \\X' U + Z' v \\ < 2 we find 



(75) 



< a. Simi- 



'« 1^') 2 

Using these facts, the triangle inequality, and (-D^) 2 = /, we obtain 



JJ> X 'u Z 'u |,/,/\ _ PC±f^t) 



1-00 



< ^2a. 



I^>- 



1^) 



ii^xi^)+^:i^)iii 



< 2(1 + v 7 ^) ^2V2~e + ^2^+ (1 + 2^)6^ < 26e2 



(76) 



with the last inequality valid for e < 1. 
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Using this estimate, and working through the estimations as in the pre- 
vious section, we find that we may set 

6 = (2/(2n + \E\) + n)yfe + 13 (n + -\E\)&. (77) 

For a simpler expression, we may use / < n and \E\ < n 2 , obtaining 

5 = (2n 3 + 4n 2 + n)^+13(^n 2 + n)e*. (78) 

Again, we may find a better estimate with more information about the 
graph. For cluster states, which have a square lattice graph, we have \E\ < 
An. We may also perform D u measurements on all vertices and set I = 0. In 
this case we may set 5 = riy/e + 39ne^. 

4 Discussion 

4.1 Estimating expected values 

The main results concern expected values, rather than experimental out- 
comes. So in order to make use of these results in any practical imple- 
mentation we must estimate the expected values using data collected from 
experimental outcomes. The obvious approach of sampling the devices many 
times and applying a Chernoff bound is problematic. In particular, we do not 
wish to assume that separate uses of a device are independent and identically 
distributed since these assumptions would be untestable and likely false in 
many practical experiments. 

One approach to this problem is that used by Pironio et al. in |PAM + 10j . 
There the authors construct a martingale, which is a sequence of random 
variables with certain properties. In particular, the random variables need 
not be independent. This allows them to use Azuma's inequality, which 
gives good bounds for martingales on how far away a sample may lie from 
the expected value without relying on independence assumptions. A similar 
approach is viable here and a preliminary analysis suggests that good bounds 
are achievable. 

4.2 Graph state computation 

Graph states are particularly interesting for their role in measurement based 
quantum computation (MBQC, [RBOlj ). In this paradigm a graph state is 
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measured, vertex by vertex, in particular bases. Each measurement may be 
interpreted as performing a unitary on a logical qubit. The composition of 
these unitaries performs a logical circuit on the logical qubits. 

A natural question to ask is whether a self-tested graph state could be 
used for MBQC to allow for self-tested computation. Unfortunately MBQC 
depends on measurements in the X-Y plane and the measurements tested 
here are all in the X-Z plane. However, the techniques used in |MM10j could 
easily be adapted to allow testing of X-Y plane measurements which would 
then allow self-tested MBQC. In fact, in the exact case the techniques used 
in |MM10] can be used with minimal changes. A preliminary analysis of 
robustness suggests that the errors scale similarly to that of lemma S] here. 
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A Classical hidden variable model for bipar- 
tite graph states with X and Z measure- 
ments 

Let G be a bipartite graph and \if)) the corresponding graph state. We give 
a local hidden variable model that is consistent will all measurements which 
are tensor products of X and Z on this state. 

We construct a local hidden variable model by randomly choosing a value 
±1 for Z' v for each v in the graph. We then set X' v to be 

K=U Z'u- (79) 

uGN v 

Now we show that this is consistent with all possible tensor product X and 
Z measurements on Let M = X Z T , Sf]T = be such a measurement. 
First, suppose that ±M can be written as a product of stabilizers of 
Using lemma [T] we have 

M = X s Z N ^ = (-1)1^)1 J] S v . (80) 

Note that, by assumption, M has only X and Z factors, so each v G S must 
have an even number of neighbours in S. Then the induced subgraph on S 
is Eulerian and we can partition the edges of the subgraph into cycles with 
no common edges (see Diestel |Diel0j for a proof). Suppose that l-E^S 1 )! is 
odd. Then there must be at least one odd cycle in this partition and then S 



20 



has an odd cycle and so does G. Since G is bipartite this must not be the 
case and in fact ^(S 1 )! is even. Hence M = Ylxes &v an d (01 M \iJj) = 1. By 
construction M' = X' S Z' N ^ = H veS X' V Z' N * = 1 and the expected value of 
M' matches that of M. 

Now suppose that M is not a product of stabilizers of \ip). Then M 
must anti-commute with at least one stabilizer and hence (ip\M\ip) = 0. 
Meanwhile, by construction 

M> = X' S Z' T = Z' N ^Z' T -. (81) 

If N(S) = T then M is in fact a product of stabilizers. This is not the case, 
so there is at least one Z' v in the above equation which is not cancelled. Since 
all the Z' v s are chosen randomly, the product of the Z' v s not cancelled will 
also be uniformly random. Thus the expected value of M' is 0. 
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